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I review recent developments in lattice calculations of B decay matrix elements and other related quantities. 



1. Introduction 

B decay phenomenology is a rapidly growing 
area of particle physics, as we expect that the 
precise measurements of various B meson decays 
provide rich information of the Standard Model 
and the physics beyond it. Following the suc- 
cessful CLEO, LEP and CDF experiments, the 
next generation B factories (BaBar and Belle) 
have just started operation, and other experimen- 
tal projects will follow. 

Lattice QCD may contribute to this program 
by calculating the B meson decay matrix ele- 
ments starting from the first principles. Its appli- 
cation is extending from the decay constant fg 
and bag parameter Bb to various semileptonic 
decay form factors. Results are being checked by 
several groups using different lattice formulations 
of heavy quark, and systematic errors are care- 
fully estimated. In addition, studies of new appli- 
cations have been started, such as the zero recoil 
form factor of B — > D^*Hu, the width difference 
of B s meson, and the lifetime ratio of b hadrons. 
In this talk I review recent developments in such 
efforts. 

This talk is organized as follows. Before dis- 
cussing the matrix elements calculation, I will 
summarize recent progress in the lattice deter- 
mination of b quark mass in the next section. 
Semileptonic decays are discussed in Section [| 
These are divided into two parts: heavy-to-heavy 
transitions (B — > D*lu, Section 3.1 
to-light decay (B — ► nlu, Section 
describe the B — B mixing, where fg and Bg 
are relevant. Recent unquenched calculations of 
fg is summarized in Section 4.1, and some up- 
dates on Bg is discussed in Section 4.2. Finally 



3.2 



and heavy- 
I then 



I will discuss new applications of lattice QCD, 
which include the calculation of B* width dif- 



ference (Section 4.3) and 6-hadron lifetime ratios 
(Section 4.4). Those calculations involve similar 
matrix elements as in Bg. A summary of the sta- 
tus as of 1998 was given by T. Draper at the last 
lattice conference 0. 

2. The b quark mass 

Although the b quark mass is not a matrix ele- 
ment related to any exclusive decay of B meson, 
it is necessary in the calculation of many inclu- 
sive decay rates. It has become clear, however, 
that there is an ambiguity in the definition of pole 
mass of heavy quark itself (renormalon ambigu- 
ity), and its determination cannot be better than 
of order Aqc d @ • In order to avoid the ambigu- 
ity one has to use a short distance mass definition, 
such as MS quark mass mf,(mf,). 

In the lattice QCD calculation, important 
progress was made by Martinelli and Sachrajda 
||, when they obtained the two- loop coefficient 
in the pcrturbative matching between continuum 
QCD and lattice HQET (static approximation) 
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where A = Mg — £ is obtained with the "bind- 
ing energy" £ measured in lattice simulations. In 
this expression the perturbative expansion in the 
first parenthesis matches the energy shift in the 



5.2 
5.0 
| 4.8 
J4.6 

-a 

S 4.4 
4.2 



"i 1 i 1 r~ 

Maitinelli-Sachiajda 



HA 



• 


2 loop 




o 


1 loop a f m b ) 




A 


1 loop a^l/a) 




V 


1 loop aJWa) 




> 


1 loop O-JUg 





y- 



0.0 0.2 



_J , I , L_ 

0.4 0.6 0.8 

a (GeV -1 ) 



1.0 



Figure 1. Lattice calculation of the b quark mass. 
Martinelli-Sachrajda's || two-loop results (filled 
circles) are plotted together with the correspond- 
ing one-loop matching results with various defini- 
tions of coupling constant (open symbols). APE's 
HQET result §1 with N F = 2 (square), and 
NRQCD result g] with N F = (star) and with 
Np = 2 (plus) are also shown. 



lattice HQET to the heavy quark pole mass. A 
power divergence exists in the loop integral that 
cancels with power divergent behavior of nonper- 
turbatively calculated A. The last line connects 
the pole mass to the MS mass. The renormalon 
ambiguity exists in both of these two expansions, 
while it cancels between them. With the pertur- 
bative expansion to two-loop, these cancellations 
are expected to become more precise, making the 
result more stable. 

Numerical results with the two-loop matching 
are very stable among three (3 values (6.0, 6.2 
and 6.4) as shown in Figure |l|, where correspond- 
ing one-loop results are also plotted with various 
definitions of the coupling constant. This result 
suggests that the cancellation of the power diver- 
gence and the renormalon ambiguity is working 
fairly well at two-loop level. 

A dynamical (Nf = 2) result along this line 
has been reported at this conference by the APE 
collaboration (V. Gimenez Q), which is shown by 
a square in Figure [l| Their result does not show 



significant shift from the quenched results. 

The NRQCD collaboration presented a b quark 
mass calculation using NRQCD action at Lattice 
98 They use a one-loop matching between 
the kinetic mass in NRQCD and the A1S mass, 
where no power divergence appears in the pertur- 
bative calculation and the one-loop coefficient is 
not large. Their results with (plus) and without 
(stars) two-flavour dynamical fermions are con- 
sistent with the HQET calculations. 

The best estimate from the lattice calculations 
rn-bi'nT'b) — 4.26±0.11 GeV, which I take from the 
two dynamical results, is in good agreement with 
recent continuum calculations 4.25(8) Q, 4.20(6) 
@, 4.20(10) g on the T resonances, and with 
DELPHI'S measurement 3.91(67) GeV at Z pole 



3. Semileptonic decays 

The lattice calculation of B meson semileptonic 
decay form factors may be used to determine \V c b\ 
and \ V u b \ with corresponding experimental results 
of exclusive decays B — > D^*Hv and B — > ■k{p)Iv. 

3.1. B -> D^lv 

3.1.1. Zero recoil form factor 

One of the most promising methods to deter- 

uwH from 
D^lv, and 



mine |V c b| is to measure \V c b\ Ti 
the differential decay rate of B 
extrapolate it to the zero recoil limit of daugh- 
ter meson w — > 1, at which the form factor 
^b^£)(»)(1) is normalized to unity in the heavy 
quark mass limit. Theoretical calculations of the 
power correction to the heavy quark limit have 
been made using QCD sum rule, but their uncer- 
tainty is still quite large (~ 9% for Tb^d-{^)) 
|Io| . Lattice calculation could be an important al- 
ternative, if it provides determination better than 
0(5%). 

Recently, the Fermilab group has proposed a 
method to calculate the deviation from the heavy 
quark limit by actually measuring the heavy 
quark mass dependence of the form factors on 
the lattice |l|,|l^]. Their key observation is that 
most of statistical and systematic errors cancel in 
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a ratio of matrix elements ll 111 



(l)l 2 - 



(D\V ch \B)(B\V bc \D) 
(D\V^\D){B\V bb \By 



(2) 



where h+~~* D (l) denotes a zero recoil form factor 
of B — ► Dlv decay through a vector current V£ b . 
The statistical error is less than one per cent for 
typical values of rrif, and m c , and thus it is possible 
to study its mass dependence. The renormaliza- 
tion factor to match the heavy-heavy current on 
the lattice to its continuum counter part is per- 
turbatively calculated, depending on the heavy 
quark mass, and found to be small for the above 
ratio Q. 

They fit their data for h + (l) as well as for /ii(l) 
and h which are form factors of B* — ► D* 
and B — > D* modes obtained through similar ra- 
tios, with the form of XjvriQ expansion predicted 
by the heavy quark symmetry Jl4| , ^5[ 
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The parameters lp, lv and A of 0(l/mg) terms 
are determined with the fit. Their preliminary re- 
sult for B — > D*lv form factor obtained at (3=5.7 
is .Fb^d-CI) = 0.935(22)(i^)(8)(20), where er- 
rors represent statistical, mass determination, 
perturbative and unknown 0(1/ttiq) in the given 
order p2| . The systematic error associated with 
the Fermilab formalism of heavy quark Jl6| is a 
subtle issue, especially because the 1 /tUq terms in 
the effective Hamiltonian are not correctly tuned 
with the use of the clover action, while the correc- 
tion to be measured is of 0(1/toq). Nevertheless, 
using the 1/toq expansion, it is generally shown 
that only 0(\/mQ) terms in the action and cur- 
rents contribute to the ratios calculated above, 
and thus the systematic error is well under control 
|l7| . The result is consistent with the recent QCD 
sum rule calculation 0.89±0.08 fllOfl, and the es- 



timated error is already slightly smaller, demon- 
strating a possibility to improve the determina- 
tion of \V c b\ using the lattice calculation in near 
future. 

3.1.2. Shape of the form factor 

In recent high statistics experiments, the ex- 
trapolation to the zero recoil limit without the- 
oretical constraints seems good enough fig)] . It 
is, however, informative to compare the shape 
of form factors with theoretical predictions in 
order to check the reliability of the extrapola- 
tion, and also to check the theoretical methods. 
Two new calculations of heavy-to-heavy form fac- 
tors (Isgur-Wise function) have been presented by 
UKQCD Jl|] and by Hein et al. §(J at this con- 
ference. 

UKQCD used the non-perturbatively O(a)- 
improved action at /3=6.0 and 6.2, and calculated 
D — * D form factor h + jl9| . They found no signif- 
icant dependence on the heavy quark mass, and 
their preliminary result for the slope parameter is 
Pl, d - 1-10(±?S)(±I) and 1.12(t*j)(tf) at /?=6.0 
and 6.2, respectively. The precision is the best 
among previous results |2l| ] , and the scaling with 
two lattice spacings is remarkable. 

Hein et al. presented a very preliminary study 
with a NRQCD action at /3=5.7 @. Form factor 
h + is obtained for B — > D and also for B — > D' 
mode, where D' denotes a radially excited state 
of D meson. Perturbative matching factors re- 
cently calculated by Boyle |^2| are incorporated 
for heavy-heavy (axial-)vector currents. 

3.2. B -> irlv 

Exclusive scmilcptonic decays B — ► Tr(p)lv 
could be used for the determination of \V u b\, pro- 
vided that corresponding form factors are theo- 
retically calculated. Unfortunately, in the lattice 
simulations, it is difficult to put large recoil mo- 
mentum on the daughter meson, so that momen- 
tum transfer squared q 2 is restricted in the large 
q 2 region (small recoil momentum). The lattice 
calculations can be, however, still useful, if statis- 
tics in the experiment is precise enough to mea- 
sure the partial decay rate for the large q 2 region, 
and such work has already been done by CLEO 
for B plv pi. 
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There are two form factors f + (q 2 ) and f°(q 2 ) 
involved in B — > irlv. f°(q 2 ) has a negligible 
contribution to the physical decay rate, as it is 
proportional to the lepton mass. It is, however, 
interesting to study f°(q 2 ), since the soft pion 
theorem relates f°(q 2 nax ) to fs, which may be 
used for a test of the lattice calculations. On 
the other hand, the calculation of f + {q 2 ) is of 
practical importance, as it is directly related to 
the physical decay rate and may be used for a 
precise determination of \V u b\ of 0(10%). 

3.2.1. f(q 2 ) 

Based on the LSZ reduction formula and cur- 
rent algebra, the soft pion theorem predicts 



2 ) 

max J 



fs 



(4) 



and p,,. — ► 0, thus q, 



2 

max 



in the limit m, 
m%. 

It was pointed out by Onogi at Lattice 97 that 

'"'"max) and Ib do not 
4|. The inconsis- 



the lattice results for / u (g 
seem to satisfy this relation 
tency has become even clearer with new data 
as shown in Figure 0, in which I plot currently 
available results for f°(q max )VM (with renormal- 
ization group correction), comparing them with 
representative data for (/b//tt) VM- The heavy 
quark scaling law predicts both quantities should 
be constant, up to 1/M corrections. 

Possible cures for this problem have been pro- 
posed at this conference: nonperturbative renor- 
malization and a method of chiral and q 2 extrap- 
olation, which I discuss in the following. 

The JLQCD collaboration has calculated a ra- 
tio of renormalization factor {Za/Zv) hl nonper- 
turbatively, using the chiral Ward identity |pl| 



Z A Z$ L I d 4 y((d^, 



2m q P){y)V HL {x)O) 



= -Z» L (A» L (x)0) 



(5) 



where and P denote the light-light axial- 
current and pseudoscalar density defined on the 
lattice, and Za represents the renormalization 
factor for A^ available nonperturbatively [34| . 
Vp ffL and Aq L are temporal component of heavy- 
light vector and axial- vector currents, which ap- 
pear in the definition of f(q%, ax ) and of fs 
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Figure 2. Check of the soft pion relation. Data 
for /°(g*lare from ELC f| UKQCD ||, 
Hiroshima JLQCD ||, Fermilab [||. Also 
plotted are new results reported at this confer- 
ence by APE @ and JLQCD jn). Light quark 
mass is fixed at the strange quark mass for some 
of the above data, since there is no significant 
light quark mass dependence found so far. Two 
representative results are taken for Jb/ f-n from 
JLQCD with clover action (open circles) |3^] and 
with NRQCD action (open squares) |33]. 



respectively. Their simulation method follows 
that of Maiani and Martinelli |55|, while they 
use 0(a)-improved action and current for light 
quark and static heavy quark action. A prelimi- 
nary result for the static-light currents at /3=6.0, 
(Za/Z v ) hl =0.72(1), is significantly smaller than 
the one-loop value 0.86(4). This result indicates 
that the one-loop calculation of the renormaliza- 
tion factor of the static-light current may contain 
a large systematic error (~ 20%), and suggests 
larger /°(gm a2; ) and/or smaller fg in the static 
limit, which gives right direction to satisfy the 
soft pion relation, but magnitude is not enough. 
An ongoing study to calculate Zj[ L nonperturba- 
tively using the Scrodinger functional method has 
also been presented at this conference by Kurth 
and Sommer pal. 



UKQCD has pointed out that a term propor- 
tional to TOjr is necessary in the chiral extrapola- 
tion as well as the usual m 2 term |37) . That de- 
pendence appears solely from the m n dependence 
otq'Lax'- Imax = ijn B -m^) 2 ~ m 2 B -2m B m„. As 
a result the chiral extrapolation could be subtle, 
since m-n and m 2 behave similarly in the region 
where simulations are made. In order to remove 
TOtt term they interpolate f°(q 2 ) to several fixed 
q 2 values for each of the active and spectator light 
quark masses. Then, the chiral extrapolation can 
be done at each fixed q 2 with the m 2 term only. 
The range of available q 2 for all light quark mass 
is substantially lower than the physical <j^ aa ., so 
they employ a pole dominance model [] to obtain 
the physical f°{qmax) from data points. This ex- 
trapolation to the physical g„ aa . makes f°(qmax) 
significantly high 1.3(i|)(-2)j which is consistent 

with / B // 7r =1.35(±g)C±^ 1 ). 

Although it is a nice observation, the method 
is, to some extent, model dependent. Therefore, 
a cross check seems necessary with the direct ex- 
trapolation method including and m 2 terms, 
which requires high statistics data at several (ac- 
tive and spectator) light quark masses. 

MILC collaboration have started such a study 
using fatlink clover quark action for both heavy 
and light quarks on their dynamical quark config- 
urations |39[] . They presented a very preliminary 
result that soft pion theorem is satisfied when ex- 
trapolated including m n term albeit with a large 
error. 

3.2.2. f+(q 2 ) 

Model independent calculation of f + (q 2 ) has 
great phenomenological importance, since the dif- 
ferential decay rate dT(B — ► nlv)/dq 2 is pro- 
portional to |U u b| 2 |/ + (g 2 )| 2 and its measurement 
may be used for the determination of \V u b\. 



' , where lattice calcula- 



tion is most powerful, the pole dominance picture 
should give a good approximation, as the corre- 
sponding pole of B* meson is very close. Using 
the B*Btt coupling g B - Br: and B* meson decay 



1 Strictly speaking they combine /+ and f° when they 
fit with dipole-pole ansatz [Eq| or with Becirevic-Kaidalov 
model fesl]. The functional form of f°(q 2 ) are pole-type 
in both models. 
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Figure 3. Heavy quark scaling of the residue of 
the pole dominance model. 



constant f B * 1 defined by 
(B(p)Tr(k)\B*(p>)) = g B * B „(k-e), 



K\B* 



= if B *m B *e, 



the form factor f + (q 2 ) is given by 



hgB'BnfB'mB' 



(6) 
(7) 



(8) 



up to the contribution of higher excited states, 
which is relatively small for large q 2 . 

Many previous lattice calculations supported 
or assumed the shape of the pole model p5| , ^o[ 
p6|]2lf |. At this conference JLQCD has further 
tested the heavy quark scaling of the pole model 
using the NRQCD action Q. They fitted ob- 
served the measured form factor with (^) and ex- 
tracted the numerator (pole residue Res/ + ). The 
heavy quark scaling law predicts gs-B-K ~ M and 
f B , ~ M -1 ' 2 , so that the pole residue should 
behave as M 3 / 2 .^ They confirmed this behavior 
as shown in Figure [|, and obtained a very pre- 
liminary result g=0.33(4) for the B*Bir coupling, 
where g is defined through g B * B7r = (2m B / 'fn)g- 

A direct lattice calculation of the B*Bir 
coupling has recently been made by UKQCD 

2 This scaling law is compatible with the scaling predicted 
for f^iq^ax) itself (~ Af 1 / 2 ), if the daughter pion mass 
is kept finite. The denominator of (^) behaves as (rrVg„ — 
m? B ) + 2mg m n ~ 2Mm, . 
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collaboration |41|j . They use the reduc- 
tion formula to relate the matrix element 
(B(p)-K(k)\B* (p 1 )) to a 'semi-leptonic transition' 
amplitude (B(p)\A fJi \B* (p + k)), where is a 
light- light axial current. Then they use the 
stochastic propagator to evaluate the latter ma- 
trix element. Their result g=0.42(4)(8) is con- 
sistent with a phenomenological determination 
through D* — ► Dir [[l2| and other phenomeno- 
logical model calculations |Q. 

In the other limit q 2 — 0, where the pion re- 
coil momentum is large, the light cone scaling 
law / + (0) ~ M -3 / 2 holds, which is not compati- 
ble with the pole dominance model that predicts 
~ A/ -1 / 2 . Becirevic and Kaidalov proposed a 
model to parametrize f + (q 2 ), introducing a term 
representing the effects of higher excited state 
contributions pq |. By choosing a parameter of 
the new term it is possible to produce an ad- 
ditional suppression of ~ M _1 , and the model 
becomes consistent with the light cone scaling 
law. Using this and other models, such as the 
pole-dipole model, one may extrapolate the lat- 
tice data to q 2 = 0. APE @ and UKQCD || 
presented / + (0) = /°(0), which are consistent 
with each other and with previous calculations 

As I discussed before, the comparison with ex- 
periment and extraction of \V u b\ is possible with 
partial decay rate without introducing model de- 
pendence. Fermilab group proposed to compare 
the differential decay rate dT/d\p n \ in the region 
400 MeV < \p n \ < 850 MeV, where systematic 
error is minimized p9| . An update was reported 
at this conference p3, where they found good 
scaling between /3=5.7 and 5.9. 

4. B — B mixing 

The mass difference of two neutral B& mesons 
is measured quite precisely AMj = 0.481±0.017 
ps _1 , and a bound is known for B s meson AM S > 
14.3 ps" 1 (H). To extract the CKM matrix el- 
ements \V t d\ and |Vt s /V t[ i|, hadronic parameters 
!%Bb and /| BbJ 'fgBs must be obtained theo- 
retically, for which the lattice calculation has been 
proven to be the best tool. Here I describe up- 
dates on the calculation of these quantities. 



0.250 



% 0.200 



JLQCD 97 
Fermilab 97 
APE 99 
UKQCD 99NP 
UKQCD 99TAD 



□ GLOK98 
O JLQCD 99 
A CP-PACS 99 



0.150 ■ 



0.10 



,Q| , 1 1 1 , 1 , L 

0.0 0.2 0.4 0.6 



1.0 1.2 1.4 l.i 



a(GeV ) 

Figure 4. Recent quenched lattice calculations 
of fg using 0(a)-improved actions. Results with 
the Fermilab formalism of heavy quarks |l6| are 
given by filled symbols: JLQCD Fermilab 
j46|. Shaded symbols represent calculation in- 
volving an extrapolation in heavy quark mass: 
APE (47|, UKQCD NRQCD results are 

given by open symbols: GLOK (49), JLQCD @, 
CP-PACS M. 



4.1. f B 

In the quenched approximation, the lattice re- 
sults have stabilized very well 0; many groups 
agree with each other using different formula- 
tion/approach for heavy quark and different lat- 
tice spacings as shown in Figure ^, where recent 
quenched results with 0(a)-improved actions are 
plotted as a function of a. The MILC collab- 
oration has also done an extensive study with 
unimproved action, and their result extrapolated 
to the continuum limit is consistent with the 



above improved action results 51 1. The dominant 
source of systematic errors (apart from quench- 
ing) depends on method and parameters used, 
but their typical size is about 10%. 

Then a natural question is what happens with 
unquenching. Before this conference, two groups 
(MILC HJ, Collins et al. ||) reported dynam- 
ical (Np = 2) simulations, and both suggested 



increase of Jb with unquenching.^] 

At this conference MILC updated their calcu- 
lation with three new dynamical lattices at ,9=5.6 
(three different sea quark masses), which is con- 
sistent with their previous dynamical result [p5[ . 
A problem in their result is a large a dependence 
seen in the quenched data due to the use of unim- 
proved Wilson quark. A linear extrapolation to 
the continuum limit gives a substantially lower 
value compared to the data at finite a. On the 
other hand, their dynamical results do not show a 
similar a dependence and the continuum limit re- 
mains high. For this reason, although their result 
suggests fg F=2 > fs F=0 , the conclusion is not 
solid enough. Therefore, they started a new cal- 
culation using the fatlink clover action for heavy 
quark, with which scaling behavior is expected to 
be improved. Q A preliminary result favors lower 
value of fs, albeit with large statistical error. 

The CP-PACS collaboration presented two new 
calculations of Jb on their dynamical lattices 
(Aj?=2)generated with an RG improved gauge 
action |5q |: 

Ali Khan discussed a NRQCD calculation at 
/3=1.95 (1/a ~ 1 GeV) with two sea quark masses 

57|. A correction of order a s /M of the heavy- 
light current induced by the operator mixing [p8| 
is included. f\ They compared the dynamical re- 
sult with their quenched result obtained at a sim- 
ilar lattice spacing, and found clear increase with 
unqucnching (~ 15-20%), while no difference was 
found between two sea quark masses. 

Shanahan presented another systematic study 
of unquenching at three (3 values (l/a=0.7~1.7 
GeV) with four sea quark masses (for each j3) 
|6l| . They used the 0(a)-improved relativistic ac- 
tion for both heavy and light quarks and applied 
the Fermilab reinterpretation for heavy quark 

16 1 . The chiral extrapolation was performed with 
m sea = m va iencei so the real unquenching was 
made for the first time. 
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Figure 5. Dynamical lattice calculations of 
Results are from MILC ^'11' Collins et al. (52 1 
and CP-PACS |57]|3l|]. Quenched results as shown 
in Figure ^ are also plotted with small open sym- 
bols. 



All dynamical results for fs are plotted in Fig- 
ure H together with the recent quenched data. We 
observe clear upward shift of fs with the inclu- 
sion of dynamical quarks. Although it is a dif- 
ficult task to combine the results from different 
groups, we can crudely say that all available data 
is consistent with the following estimates: 



!b (MeV) 
fB. (MeV) 



N F =2 N F =0 

210 ± 30 170 ± 20 

245 ± 30 195 ± 20 

1.16 ±4 1.15 ± 4 



where I also list the results for fs B and fs, / fs- I 
do not attempt to extrapolate these results to the 
physical Np = 3 limit. To do so, it seems neces- 
sary to understand the systematic errors coming 
from the use of different actions and lattice spac- 
ings. The sea quark mass dependence should also 
be clarified. 



3 Other theoretical studies, such as Quenched Chiral Per- 
turbation Theory Jjjl , quark model, 'bermions' jjjj ], sug- 
gested the same conclusion. 

4 Necessary one-loop calculation has been performed by 
Bernard and DeGrand [159]. 

5 Necessary one-loop calculation with the RG improved 
gauge action has been performed by Ishikawa fcol. 



4.2. B B 

In contrast to the achievement for fs, the lat- 
tice calculation of Bb is still premature. 

In the static approximation, the 0(a)-improved 
results by Gimenez and Martinelli |32| and by 
UKQCD [p3| have been reanalyzed in a recent 
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paper by Gimenez and Reyes |64j using corrected 
one-loop matching coefficient, and a disagree- 
ment, which existed between the static-clover 
and static- Wilson results |66fl, has been greatly re- 
duced]] 

The Hiroshima group performed a calculation 
using the NRQCD action I^Q, and found sig- 
nificant decrease of B B (m,b) as one includes 1/M 
corrections. The matching, however, was done 
with the coefficient in the static limit, and thus 
large 0(a/(aM)) systematic error is expected for 
the slope in 1/M. 

In the calculation with relativistic actions, 
UKQCD presented the first calculation with the 
0(a)-improved action at /3=6.0 and 6.2 at Lattice 
98 |H , which has recently been updated Q . To 
obtain the result at the B meson mass, an extrap- 
olation from charm mass regime is necessary, and 
they found a clear negative slope in 1/M. 

At this conference, APE group Q has pre- 
sented the first result obtained using nonpertur- 
bative renormalization |7l|| . They found a similar 
dependence of on 1/M, but their final nu- 
merical results are not yet available at the time I 
wrote this contribution. 

Figure ^| presents a compilation of lattice data 
for <P BB (pb) = (a s (Mp)/a s (MB)) 2/l3o B B (p b ) 
with pb—5 GeV as a function of 1/Mp. The 
renormalization factor is introduced to cancel the 
ln(M/fib) dependence appearing in the match- 
ing factor fig . It is encouraging that all rela- 
tivistic results including the early works 72 7j| 
show a reasonable agreement with each other, 
and that the recent UKQCD data show a nice 
scaling between /3=6.0 and 6.2. The extrapola- 
tion to the static limit (~ 0.92), however, seems 
considerably higher than the 0(a)-improved re- 
sults in that limit. It suggests that there is un- 
known sources of systematic error in either or 
both of static (NRQCD) and relativistic calcula- 
tions. Higher order perturbative corrections (for 
both) and 0((aM) 2 ) uncertainty in the relativis- 
tic calculations are their potential candidates. 
For this reason, my summary of the current avail- 



6 UKQCD has also presented a corrected analysis 
But I do not use their number in this summary, since TEe 
tadpole improvement is not performed in their analysis 
and there seems large systematic error remaining. 
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Figure 6. 1/M dependence of <& Bb (5GeV). The 
static and NRQCD data are from Kentucky j(K]], 
Gimenez-Reyes @ (Reanalysis of §§]), UKQCD 
|H (Reanalyzed by Q), and Hiroshima 
The relativistic calculations are Bernard et al. 
\7%, ELC f7|, Bernard-Soni [Q, Gupta et al. 
[7|, and UKQCD §|. Open symbols are ob- 
tained with Wilson quark for heavy and/or light 
quarks, and filled ones are 0(a)-improved. 



able data includes a large systematic uncertainty: 
B B {m b ) = 0.80(15). 

The allowed region on the (p, r/) plane of the 
CKM matrix is shown in Figure |7|. The two 
flavour result f B = 210 ± 30 MeV and the conser- 
vative estimate of B B are used to draw the con- 
straint from AMrf. Due to the upward shift of f B 
from the previous quenched results, the allowed 
region favors p > 0. 



4.3. B s width difference 

The width difference in the B s 
given as 



B s mixing is 



Ar s cx 

Im- 



2M Bs 



B s \i / d A xTH ca {x)H cS mB s )i$) 



where H cS represents the AB—1 effective Hamil- 
tonian |76| . Only the final states into which both 
of B s and B s can decay contribute. The l/mt ex- 
pansion induces two four-quark AB=2 operators, 
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Figure 7. Constraint on the CKM matrix ele- 
ments. 



whose matrix elements with B s and B s states are 
Bb and Bs, where Bs is dehned through 



(B a \O s (ji)\B. 
5 



-J. 



2 Ml 



(nib + m s ) z 



(10) 



and O s = b(l - 7 5 )s6(l - 7 5 )s. 

At this conference, the Hiroshima group pre- 
sented a calculation of Bs using the NRQCD ac- 
tion Q . Their calculation method is the same as 
that of B B and they obtain S s (m 6 )=1.19(2)(20). 
Using a next-to-leading order formula of Bcnckc 
et al. W%, the width difference is obtained as 
(Ar/r) b s = 0.16(3) (4), where errors are from fg s 
and Bs respectively. The two-flavour result for 
fs, discussed in Section 4.1 is used. The latest ex- 
perimental bound from DELPHI is (Ar/r) B , < 
0.42 |78 . 



4.4. Lifetime ratios 

The ratios of lifetime of b hadrons, such as 
t(B-)/t(B°) and r(A b )/r(B ), provide an im- 
portant test of the theoretical method to calcu- 
late the inclusive hadronic decay rates [jr9[ . In 
the 1/mt, expansion, the leading contribution to 
the decay rate comes from a diagram in which 
the &-quark decay proceeds without touching the 
spectator quark, so that it does not contribute to 



the lifetime ratios. The 0(1/ ml) correction to 
the ratios is also small for the same reason, and 
the first correction involving the spectator quark 
effect is of 0(1/ ml), which is parametrized by 
the 'B parameters' of AB—0 four-quark opera- 
tors. UKQCD computed these matrix elements 
for the first time and obtained t(B~)/t(B°) = 
1.03(2)(3) ^35|, which is consistent with the re- 
cent experimental result 1.07(2) po| . 

It is a known problem that the lifetime of A& 
is surprisingly shorter than that of B mesons 
T(A b )/T(B°) = 0.79(5) @. It is, therefore, in- 
teresting to see whether it is explained with the 
theoretical calculation, in which the similar ma- 
trix elements of four quark operators for Af, are 
required. The UKQCD group has studied these 
matrix elements at (3=5.7 with 20 gauge configu- 
rations |8l| , and found that the spectator effect is 

large 6%. Although their result r(A 6 )/r(S ) 

= 0.91(1)^0.93(1), depending on the light quark 
mass, is much higher than the experimental value, 
higher statistics calculations at higher (3 values 
seem necessary to draw a definite conclusion. 

5. Conclusions 

Lattice calculations provide model independent 
predictions for many important B decay ma- 
trix elements. Progress made for the zero re- 
coil B — > D^lv from factors is essential for pre- 
cise determination of \V c b\, and the shape of the 
form factors is also being studied with improved 
techniques (NP improved action, NRQCD, etc.). 
More study is necessary to achieve a complete 
understanding of the B — ► nlv form factors: the 
violation of soft pion theorem and the shape of 
f + (q 2 )- The determination of \V u b\ with 10% pre- 
cision will become possible, once we understand 
these questions. 

The dynamical quark simulation has become 
practical by several groups, and its effect on fs 
has been identified. Further systematic study 
like that of MILC and CP-PACS is necessary to 
understand systematic errors and eventually to 
obtain physical prediction at Np = 3. An un- 
quenched study of other quantities should also be 
important. 

Several new applications have also been stud- 
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ied, such as the width difference of B s meson, 
and the lifetime ratios of b hadrons. Those 
will become useful theoretical calculations, pro- 
vided that statistical and systematic errors are 
improved. 
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